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We explore if the phase separation or coexistence of the spin density wave (SDW) and supercon- 
ductivity (SC) states has any relation to the incommensurability of the SDW in the Fe pnictide 
superconductors. A systematic method of determining the phase separation or coexistence was 
employed by computing the anisotropy coefficient /? from the the 4th order terms of the Ginzburg- 
Landau (GL) expansion of the free energy close to the tricritical/tetracritical point. It was comple- 
mented by the self-consistent numerical iterations of the gap equations to map out the boundaries 
between the phase separation and coexistence of the SDW and SC phases, and between commen- 
surate (C) and incommensurate (IC) SDW in the temperature-doping plane. Our principal results 
for the sign reversed s-wave pairing SC, in terms of the multicritical temperature, Tc, the phase 
separation/coexistence boundary between the SDW and SC, T* , and the boundary between C/IC 
SDW, T^j, are: (a) IC-SDW and SC coexist for Tc < T* and phase separate otherwise, (b) SDW 
takes the C form for Tc > T^i and IC form for Tc < T^f, and (c) the thermodynamic first order phase 
transition intervenes in between the C-SDW and IC-SDW boundary for large T'Ij, where T'Ij is the 
SDW transition temperature at zero doping, T* = 0.35 T'Ij and T^.j = 0.56 T'Ij. The intervention 
makes the phase diagram more complicated than previously reported. By contrast no coexistence 
was found for the equal sign pairing SC. These results will be compared with the experimental 
reports in the Fe pnictide superconductors. 



PACS numbers: 74.70.Xa, 74.25.Dw, 74.25.Ha. 



I. INTRODUCTION 

The alluring prospect of opening a key window to un- 
derstanding the mechanism of high temperature super- 
conductivity has attracted enormous research activities 
in the iron based pnictidesji"— A widespread school of 
thought regarding the pairing interaction maintains that 
magnetic fluctuations are intimately involved for the su- 
perconductivity. This view seems natural from the over- 
all phase diagram of the pnictides in the temperature 
and doping plane. The superconductivity (SC) emerges 
out of the parent antiferromagnetic (AF) state as the 
static AF order is suppressed and charge carriers are in- 
troduced through the electron or hole doping in common 
with the cuprate high temperature superconductors.* 
This view is further supported by the neutron scatter- 
ing investigations showing a resonance at the wave vec- 
tor related to the AF order coincident with the onset of 
superconductivity)^ 

The intimacy of SC and AF or spin density wave 
(SDW) states is a common feature among the uncon- 
ventional superconductors like the cuprates, pnictides, 
and heavy fermion superconductors.® It will, therefore, 
be important to understand the interplay between the 
SC and SDW orders. For the cuprates, it is perhaps one 
of the best established observations that the SC emerges 
out of AF parent state as the doping is introduced.^ For 
the heavy fermion Cell 5 compounds, the interplay be- 
tween them has been studied actively for CeRhlus and 
CeCoIns. For CeRhlus, for instance, the AF and SC 
phases coexist under pressure at zero magnetic field. ^'^ 



Within the coexisting dome, the AF changes from in- 
commensurate (IC) to commensurate (C) SDW as the 
pressure is increased.^9 

Experimental investigations on the interplay between 
the SC and SDW for the pnictides report diverse results 
with regard to the pnictide families and the dopants. 
For the 1111 family, Luetkens et alt^ reported from 
the muon spin relaxation (/iSR) and Mossbauer spec- 
troscopy on the LaFeAsF^jOi-a; (Lallll) compounds 
that the magnetic state disappears and superconductiv- 
ity emerges abruptly as the doping x is increased. For 
Smllll system. Drew et al. found that the AF and SC 
regions coexist which could be due to phase separations^ 
The neutron scattering measurements for the Cellll 
compounds reported a magnetic phase diagram similar 
to Lallllji^ For the 122 family, Laplace et ali^ re- 
ported that the IC-SDW coexists with the SC on the 
atomic scale in Ba(Fei_j;Co2;)2As2 compounds by mea- 
suring ^^As nuclear magnetic resonance (NMR) and sus- 
ceptibility. Julien et al^ contrasted K and Co doped 122 
compounds by also measuring ^^As NMR and found that 
the SC and SDW phase separate for Bao.6Ko.4Fe2As2 but 
microscopically coexist for Ba(Fei_a;Co2;)2As2 in accord 
with Laplace et al. Parker et ali^ reported from the neu- 
tron and muon experiments on NaFei-ajM^^As (M=Co, 
Ni) that the SC and SDW coexist on atomic level. 

Theoretically the interplay between SC and SDW for 
the Fe pnictide superconductors was studied in a sim- 
ple two band model by many groups. ^^^^^ They reported 
that the coexistence of SC and SDW is possible when 
SDW is IC for the tt phase shifted pairing. In this pa- 



per we also take the two band model to study the inter- 
play between the SC and SDW. We make the Ginzburg- 
Landau (GL) expansion of the free energy close to the 
multicritical point for IC as well as C-SDW states. The 
nature of transitions was then determined from the 4th 
order terms, that is, whether they are continuous or dis- 
continuous and whether the two orders coexist or phase 
separate. The phase separation or coexistence was de- 
termined by computing the anisotropy coefficient P from 
the 4th order terms of the GL expansion for general IC 
SDW states. The superconducting state was modeled in 
terms of the sign reversed two band SC theoryi^"— The 
< /3 < 1 or not means phase coexistence or separation. 
See Eqs. (|35l) and ([36|) below. The multicritical tem- 
perature Tc which equals to the superconducting critical 
temperature in the simple model we took is set by the 
interaction A5. The SDW transition temperature at zero 
doping T^j = Tm{5 — 0) is set by \m- S is the deviation 
from the perfect nesting and may be tuned by doping or 
pressure. 

When Am is only slightly larger than As, then 
T^j/Tc ^ 1, and the multicritical point occurs at a small 
doping. The deviation from the perfect nesting is then 
small and the SDW takes the commensurate form. The 
computation of /3 yields /? > 1 which means that the 
C-SDW and SC phase separate. As the ratio of T^j/T, 



gap equations from a simple two band model. ^^ Then we 
introduce the anisotropy coefhcient (3 from the 4th order 
terms of the GL free energy which determines the phase 
separation or coexistence between the SDW and SC. In 
section III, we present the detailed calculations of the 
multicritical point, the coefficient (3, and SDW and SC 
order parameters. We first show the multicritical point 
in the T — 6 plane in Fig. [1] from which one can read off 
the nature of phase transitions. Then we present three 
typical cases: phase separated C-SDW and SC in Fig. 2, 
discontinuous SDW and SC phases in Fig. 3, and coex- 
isting IC-SDW and SC phases in Fig. 4. These results 
will be commented in comparison with the experimental 
situations. Section IV is for summary and concluding 
remarks. 



II. FORMALISM 

We take the following hamiltonian as with Ref.^^ to 
describe the interplay between superconductivity and an- 
tiferromagnetism. 



H=J2 ^»fecL<,c,fc<. +Hs + Hm, (1) 



becomes IC for Tc < Tlj . We use the notations that T^ 
represents the boundary between C-SDW and IC-SDW 
determined by Eq. (P^ and T* represents the phase sep- 
aration/coexistence boundary between SDW and SC de- 
termined by Eq. (I36|) below. T^j > T* in our model. The 
IC-SDW and SC phases remain separated for Tc > T*. 
For Tc<T*, the ^ becomes < /3 < 1 and the IC-SDW 
and SC phases coexist. We obtain T^j = 0.56 T"^ in 
agreement with RefJI, and T* = 0.35 T^^. See Eqs. ^ 
and (gl]) below. 

The SDW transition as T is reduced is continuous for 
almost all of the parameter space. But, as the doping at 
the multicritical point is increased, the 4th order term of 
the SDW order /?m of Eqs. (|23)) and (|A5|) becomes neg- 
ative and the first order transition intervenes in between 
the C-SDW and IC-SDW boundary. The intervention of 
the discontinuous transition makes the phase boundary 
more complicated. See the Fig. 1 and discussions be- 
low for details. Note that this 1st order transition is as 
the temperature is lowered as presented in Fig. 3. On 
the other hand, the 1st order transition between SDW 
and SC as 5 is varied below Tc corresponds to the phase 
separation as shown in Fig. 2. 

The GL expansion was combined with the self- 
consistent numerical iterations of the gap equations to 
map out the boundary between the phase separation and 
coexistence of SDW and SC in the plane of T and 5 for 
the pnictide superconductors. See the the figures 2-4. 

This paper is organized as follows: In the following sec- 
tion we will present the functional integral formulation of 
the Ginzburg-Landau free energy and the self-consistent 



increases, Tm{5) equals to Tc at a larger 5, and the SDW -H^s = 2_. ^S' [ Ci,fc,t4 -fc,iC2,-fc',iC2,fc',t + ^-C- , (2) 

becomes IC for T„ < T.*,. We use the notations that T.*, fc,fc' 

Hm = -VM^(yc\k^^C2M+q,a ^ Cr'c^,_^g ,^, Ci,fc.,^- , (3) 



where the Vm and Vs are the magnetic and supercon- 
ducting interaction strengths, respectively, and q is the 
incommensurability of SDW. q = Q and / represent 
the C and IC-SDW which is determined by maximizing 
the magnetic susceptibility x(9)- See Eq. ([201) below. 

The partition function of the Hamiltonian of Eq. ([T]) 
in the form of functional integral is given by2S (r is an 
imaginary time, and j3 = l/ftsT) 



£)[c|,4,Cl,C2 



] g- r dr[^y,drC^+cldrC2)+H[r)] ^ 



(4) 



The next step is to employ the Hubbard-Stratonovich 
transformations to decouple the electron-electron inter- 
actions of Hs and Hm- It is straightforward to obtain — 



■j^drHs(r) 



yi5[A^A;,Al,A2]e-^^ 



-Jl^drHMiT) ^ 



D[m* , m\e 



-Sa 



(5) 



where 

Ss = 



/^^? 



A^(r)A2(T) 
Vs 



{cU.tM4,-fc.,WA2(r) 



s 



M 



f J X^ \ m*{T)m{T) 

\ m(r) Y^ fJC^ ,,(r)c2,fc+g,<T(T) + h.c. \ 



Ai/2{t) are nothing but the (fluctuating) superconduct- 
ing order parameters. Ai and A2 are taken as the op- 
posite sign in accord with the s± pairingj22i2^ It is more 
convenient to introduce 



A = V-A1A2 (8) 

as with Ref<2^. Integrate out the fermions to get 



F 



dr 



d^— A*(t)A(t) 



--—m* {T)m{T) - hidet(Ms + Mm) 
Vm 



(9) 



We then make the saddle point approximation which 
is determined by the condition that the first order func- 
tional derivative of the action F with respect to A and 
m vanish. 



5F^ 



Ao,mo 



5F 
6m 



(10) 



Ao,mo 



where Aq, ttiq denote the saddle point values. These func- 
tional derivatives can be computed using the following 
matrix identity 

JlndetM = JTrlnM = Tr{M-^SM). (11) 

They are, of course, the usual gap equations:— 



--j:. 



<j=± 



diXs tanh ( — - 



(12) 






d^X 



aS 



M 



2Ea 



y^^ + TO^ 



tanh 



2T 



where the subscript was dropped, 

\s = NfVs, \m=NfVm, 



(13) 



(14) 



where Np is the density of states at the Fermi level, and 
the energies are given by 



E± 



y/m^+e±6) +A^. 



(15) 



((J^e took ^ifc — —£_ + S and ^2fc = C + ^ for the hole and 
electron Fermi surfaces, respectively.— The free energy is 
given as follows: 



F = 
(7) i- A* A + -^m*m - 2T VV In cosh (i/3^^fc 



As 



(16) 



We take Wc as the unit of energy {ujc = 1). 
The second order terms are 



F^'^ = (Am - Xm) m*m + (A^^ - xs) A*A, (17) 
where the pairing and magnetic susceptibilities are given 

by 



1 /■"'^ 1 / £ 

Xs(T) = -|^^dC-tanh(X 

XM('5,T) = iy'%£itanh(^^ 



(18) 
(19) 



When the doping is increased the magnetic transition 
may be incommensurate, that is, xm can be maximum 
for IC wave vector. For IC-SDW {q ^ 0), the magnetic 
susceptibility is given by 



XM{5,T,q) = - 



d^ 



1 



£, — q cos 9 



tanh 



2T 



(?0) 



where ( )g implies the angular average. Other technical 
details for IC case are collected in the appendix. The 
Xm(^, T, q) needs to be maximized with respect to the q 
for given 8 and T. The magnetic transition temperature 
Tm and superconducting transition temperature T^ are 
determined by vanishing second order coefficients. 

For comparison we also considered the equal sign s- 
wave pairing to check how the interplay differs from the 
sign reversed s-wave pairing. It is simple to see that the 
formula remain the same except the energies E± of Eq. 
(|15p . For equal sign pairing, the energies are given by 



E± 



£2 + A2 + J2 ± 2yJ{'m?+i'^)5'^ + m?Al2l) 



The expressions for the xSi Xm, and the 4th order co- 
efficients remain the same, except the /3' of cross term 
of Eqs. (|23p and (|A6|1 . It will alter the anisotropy co- 
efficient /3 of Eqs. JSSl) and §^ and the SC and SDW 
phases always separate for the equal sign pairing. 

Fig. 1 shows the magnetic transition temperature Tm 
as a function of S determined by 



XM{5,TM,q) = 



M 



(22) 



for Am = 0.28, 0.4, 0.52, 0.68, and 0.8. Also shown 
in the figure are T^ and T* which are the boundaries 
between C and IC-SDW, and between phase separa- 
tion and coexistence. T* can be determined from the 



0.4 



0.3 - 



t^ 0.2 



0.1 



0.0 



M 



0.68 

" 0.52 

0.40 
0.28^ 



A\ Bo . 
CV V 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 

S 



FIG. 1: The plot oiTM{5) for Am = 0.28, 0.4, 0.52, 0.68, and 



0.80, and T^ and T* . The unit of energy is ujc 



= !)• n, 



is the boundary between g = (C) and g 7^ (IC) SDW from 
solving Eq. H22p . T* is the boundary between the phase sep- 
aration and coexistence between SDW and SC determined by 
Eq. H36[). For a given Am, in terms of the multicritical temper- 
ature Tc, the region of T%i > Tc > Tli (black) corresponds to 
the phase separation between C-SDW and SC, T^ > Tc > T* 
(cyan) to the phase separation between IC-SDW and SC, and 
Tc < T* (pink) corresponds to the coexistence between IC- 
SDW and SC. As Xm is increased the region of the discontin- 
uous thermodynamic transitions intervenes around T^j and in 
the low temperature region as indicated by the green and red 
lines. The points marked by A, B, and C correspond to the 
representative cases of the phase separated C-SDW and SC, 
phase separated discontinuous C-SDW and SC, and coexist- 
ing IC-SDW and SC phases, respetively. The A, B, C cases 
are presented in more detail in Figs. 2, 3, and 4, respectively. 



anisotropy coefficient /3. It would have been extremely 
difficult to map out the complete C/IC SDW phase sep- 
aration/coexistence with the SC phase without the sys- 
tematic approach using the coefficient /3. It was com- 
puted for general q j^ IC-SDW as follows. 

The nature of the transitions, that is, whether the mag- 
netic and superconducting transitions are continuous or 
discontinuous, and whether A and m coexist or not is 
determined by the 4th order terms. Expansion of the GL 
free energy of Eq. p^ up to the 4th order yields 



F(4) =/3jv/|m|4 + /?A|A|4 + /?'|mnAp 



(23) 



All other combination terms vanish. The expressions of 
the coefficients are collected in the appendix. If all coef- 
ficients /3a/ 7 I^A, /3' are positive, then the magnetic and 
superconducting transitions as T is reduced are continu- 
ous. Depending on the relations among the coefficients 
as will be discussed below, m and A may coexist or re- 
pel each other as the doping is varied below Tc- The 
negative /3m, on the other hand, means discontinuous 
magnetic transition as T is reduced. This occurs in the 



vicinity of the C-IC transition in the present model when 
Am is large, as will be discussed below. 

In order to determine the phase coexistence/separation 
we follow the standard statistical mechanics procedure 
and focus around the multicritical point where Tm{5) = 
T^^L2» We write 



X5(T)-A-i 
Xm(T) - XjJ = Am 
where < Am < 1 is given by 



T 


-T, 




Tc 


T 


-Tc 







Am = T, 


1/ 


nUJc 


1 ^ + S 


2\ 


l./^e^ 


-qcos0 2Tc "^ 



Tc 



dxM{Sc,Tc,qc) 
dT 

2Tc 



(24) 
(25) 



(26) 



We now introduce the field r] to make the 2nd order 
terms isotropic in the order parameter space. 



(27) 



and write 



^(2) 



T-Tc 

Tr 



{H 



(28) 



Collecting the 2nd and 4th order terms we have 

A|2)+aH4 + 6|Ar + chnA|^ 



F=L^i\^\^ 



Tc 



where the coefficients are given by 
a =-2-, b = PA, c = 

The isotropic 2nd order terms suggest writing 
|77|=rcos6', |A|=rsin6'. 
The free energy can now be written as 



Am' 



(29) 



(30) 



(31) 



F = 



T -T, 



Tc 



c 2 I 4 

r + r 



[a cos* 9 + h sin"* 6 + c cos^ 9 sin^ 



(32) 



For 



a + 6-oO, (33) 

the free energy takes the minimum with respect to 9 when 

cos^ 9 = P, (34) 



where /? is given by 



/3 = 



b-c/2 
a + b ~ c 



(35) 



We then have 



m 7^ 0, A = 0, for /3 > 1, 
m = 0, A =^ 0, for /3 < 0, 
m ^ 0, A ^ 0, for < /3 < 1. 



(36) 



Eq. ([M)) . of course, represents the coexisting region. It 
may be rewritten in terms of a, 6, c as 



a > c/2 and 6 > c/2. 



(37) 



Otherwise, the two orders phase separate and the first 
order phase transition shows up as S changes below Tc- 
The coefhcients a, b, and c can be calculated from Eqs. 
(P5)) and ([50]) . Notice that the coexistence condition of 
Eq. (j37p is more restrictive than the condition 



ab > c^/4 



(38) 



2 Tc-T a- 


-c/2 


2Tc ab - 
1 Tc-T b- 


-c/2 


AIj 2T, ab- 


-c2/4 



given by Vorontsov et ali^ The value of j3 determines 
the ratio of the two condensates by Eq. (pT]) to give 



(39) 



The GL expansion is only valid close to the multicrit- 
ical point. As T is lowered further below Tc, the order 
parameters become larger and we need next-order terms 
in the GL expansion to describe the critical phenomena 
with the same accuracy. This is very tedious and in- 
efficient. We therefore obtain the order parameters m 
and A self-consistently with numerical iterations as ex- 
plained before. In the following section we will present 
the results. 
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FIG. 2: (a) The phase diagram of the SDW and SC for \m = 
0.68 and As = 0.55 where the tricritical point is at the point 
A on the second line from the top in the Fig. [l] It shows the 
phase separation between C-SDW and SC. (b) The SC order 
parameter A and SDW order parameter in as a function of 
temperature a.t 5 = 5c — 5step as indicated by the vertical 



dashed line in (a). 
function of 5 at T = 
line in (a). 



(c) The order parameter A and m as a 
0.1 as indicated by the horizontal dashed 



III. RESULTS 

To map out the phase boundary between the SDW and 
SC phases we first calculate Tm{5) which is determined 
by \m via Eq. ([22)1 . We considered the cases \m = 0.28, 
0.4, 0.52, 0.68, and 0.80 for representative values and 
corresponding Tm{5) are shown in Fig. [T] The incom- 
mensurability q was computed by maximizing xm of Eq 
(|22|) with respect to q for given T and 5. q = and 
y^ mean, respectively, the C-SDW and IC-SDW. The 
boundary between q = and g 7^ are marked on the 
Tm{5) curves by Tlj. We found 



Tl, = 0.56 T" 



M 



in agreement with Ref J^. This corresponds to 

0.958 



2r, 



(40) 



(41) 



which can be calculated from xm of Eq. ((20l) for q = 0. 

We also computed the anisotropy coefficient (3 along 
Tm{5) as the multicritical point is lowered along the 
curve. For a given Xm, a smaller A5 means that the 
multicritical point where Tm{6c) = Tc is at larger 6c and 
smaller Tc- The computation of (3 yields that there ex- 
ists the critical T* such that m and A phase separate for 
Tc > T* and coexist for Tc <T*. We found 



rji^ 



0.35 Tl,. 



(42) 



M 



In Fig. [I] the black color of each Tm{S) curve represents 
the phase separated C-SDW and SC where r°^ > Tc > 
Til, the cyan the phase separated IC-SDW and SC where 
T^j > Tc > T* , and the pink represents the coexisting 
IC-SDW and SC phases where T* > Tc. The coexisting 
SDW and SC are possible only for IC cases for the sign 
changed pairing in agreement with Refi^^. 



An interesting observation is made for large Am such 
as A = 0.8 of Fig. [TJ The coefScient Pm of Eq. d^S]) of 
the 4th order SDW term becomes negative in the vicinity 
of the boundary between the C and IC-SDW. It means 
that the SDW transition as T is reduced is discontinu- 
ous. See the discontinuous changes of the order parame- 
ters as the temperature is reduced as shown in Fig. 3(b). 
IncidentaUy, the thermodynamic first order SDW tran- 
sition was not reported in the weak couphng calculation 
of Vorontsov et alXL To check the discontinuous C-SDW 
transition as the mark B of Fig. [1] stands for, put q — Q 
in Pm given by Eq. (|A5|) . Take the limit ujc ^>- oo and 
put X = i/2T and y = 5/2T. 



A 



1 



M 



{2T)'- 



f°° dx 
■ I — [tanh {x + y) — X sech^ {x + y)](.43) 



The integral changes sign at y = 0.955 which means that 
for 5I2Tm > 0.955 the first order SDW transition oc- 
curs. Comparing this with the S/2Tlj- = 0.958 of Eq. (|iT|) 
means that the first order SDW transition intervenes near 
the C-SDW and IC-SDW boundary. The self-consistent 
numerical calculations indeed confirm this result as pre- 
sented in Fig. 3. In our calculation it does not show up 
when Am is small but begins to emerge when Xi\i > 0.5 
probably because the discontinuity for small Am is too 
small. As S increases further, however, the nonzero q in- 
creases the xm , and IC second order transition preempts 
the first order transition as shown in Fig. 1. 

Then we consider three cases where the points marked 
by A, B, and C on Fig. 1 are the multicritical point and 
calculated the phase diagram in the T — 5 plane. The 
SDW and SC order parameters m and A as a function 
of T and 5 are calculated by solving the gap equations 
via numerical iterations. The results like the phase sep- 
aration/coexistence are fully consistent with the GL ex- 
pansion. Let us first consider the point A. The compu- 
tation of /3 at the A point yields /3 > 1 corresponding 
to the phase separation between C-SDW and SC. The 
phase transition line below Tc was calculated with the 
numerical iterations of the gap equations of Eqs. (fT2|) 
and (I13p . The results are shown in the Fig. [51 Fig. 
[5](b) and (c) present, respectively, the order parameters 



A and m as functions of T at (5 



5^ 



^ step 



and as 



functions of (5 at T = 0.01 as indicated by the dashed 
lines, m shows the second order phase transition as a 
function of T, and the A and m show the first order 
phase transition as functions of S below Tc- This indi- 
cates that the SDW and SC orders phase separate as 
the doping is varied as was mentioned in the Introduc- 
tion for the 1111 compounds. Lallll,^^ Snillll,^^ and 
Cellll compounds^'' all showed the phase separation as 
the doping was changed. 

Now, let us turn to the case where the multicritical 
point is at the point B corresponding to Xm = 0.80 and 
As = 0.55. As alluded previously, the 4th order term of 
the GL expansion /3m at the point B becomes negative 
in this case and the SDW transition as T is reduced is 
discontinuous. The phase diagram in the T — S plane 
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FIG. 3: (a) The phase diagram of the SDW and SC for 
Xm ~ 0.8 and As = 0.55 corresponding to the point B in 
the Fig. [l] It exhibits the discontinuous change of transi- 
tion temperatures between C-SDW and SC at the critical Sc- 
(b) The A and m as a function of temperature at 5 = (5c as 
indicated by the vertical dotted line in (a), (c) The order pa- 
rameter A and m, as a function of 5 at T = 0.1 as indicated 
by the horizontal dotted line in (a). 



obtained by the self-consistent numerical calculations of 
the gap equation is shown in Fig. 3(a). The temperature 
dependence of the order parameters m and A is shown 
in (b). The m exhibits the abrupt onset at T = Tm- The 
experiments were reported to be the continuous transi- 
tions for 771 as a function of T.^^ The transitions, however, 
could be equally well described as a weakly first order as 
shown in the present calculations. On the other hand, 
the abrupt change of Tm and Tc as a function of doping 
was not seen in the present model calculation. The ex- 
periments reported that the SDW (SC) showed up in the 
orthorhombic (tetragonal) phase which indicates the po- 
tential importance of the spin-lattice coupling. The dis- 
crepancy between the experiments and the present cal- 
culation could be due to the neglect of the spin-lattice 
coupling. 

Fig. 4 shows the case where the tetracritical point is at 
the point C corresponding to Xm = 0.52 and As = 0.30. 



< /? < 1 is this case which means the IC-SDW and SC 
phases coexist microscopically, that is, both m and A or- 
der parameters are nonzero at the same point in the real 
space. The phase diagram of Fig. 4(a) was obtained with 
the numerical iterations of the gap equations of (jAip . 
(|A2|) . and (|A3p . The order parameters change contin- 
uously and the all the lines represent the second order 
transition. The green shaded region of Fig. 2] represents 
the parameter space where SDW and SC coexist micro- 
scopically. Figure (b) plots m and A as functions of T 
at J = 0.18 indicated by the dotted vertical line in figure 
(a). Note that m decreases as T is lowered below the co- 
existing region because SC begins to partially gap out the 
Fermi surface. Figure (c) shows m and A as functions of 
(5 at T = 0.019 indicated by the dotted horizontal line in 
the plot (a). The SDW and SC order coexist in the region 
of 0.16 ^ 6 ^ 0.21 in the low temperature limit. This re- 
sult nicely corresponds to the Co doped 122 compounds. 
As was mentioned in the Introduction, Laplace et alr^ 
reported that the IC-SDW coexists with the SC on the 
atomic scale in Ba(Fei_a;Coa;)2As2. This result was later 
confirmed by Julien et al on Ba(Fei_a:Coa,)2As2^ and by 
Parker et al. on NaFei^o^M^jAs (M=Co, Ni) J^ 



IV. CONCLUSIONS 

In this paper, we examined the interplay between the 
spin density wave and tt phase shifted superconductivity 
in the Fe pnictide superconductors. We have obtained 
the phase diagram in the plane of the temperature T 
and chemical potential 5 with the combination of the 
Ginzburg-Landau expansion of the free energy near the 
multicritical point and the self-consistent numerical iter- 
ations of the gap equations. By calculating the multicriti- 
cal temperature Tc as a function of the chemical potential 
5 as shown in Fig. 1, we presented possible cases of phase 
separation/coexistence among the commensurate SDW, 
incommensurate SDW, and superconducting phases. 

Then three typical cases were considered in more de- 
tail. The phase separation of C-SDW and SC for Tc > 
T^ was shown in Fig. 2. The phase diagram in the T ~ S 
plane, the T and S dependences of the SDW and SC or- 
der parameters were shown. In Fig. 3, the discontinuous 
SDW-SC transition case was presented. And in Fig. 4, 
the coexisting IC-SDW and SC for Tc < T* case was 
presented. 

In doing so, we employed a systematic way of deter- 
mining the phase separation or coexistence between SDW 
and SC orders from the 4th order terms of the Ginzburg- 
Landau free energy expansion. We found that if Tc is 
larger than T* then the two orders phase separate, but 
coexist for Tc < T* unless the first order transition inter- 
venes. Although the SDW transitions as temperature is 
lowered are continuous for most of the parameter space, 
they can become first order if the doping at the multicrit- 
ical point becomes large. Then the first order transition 
intervenes in between the C-SDW and IC-SDW bound- 
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FIG. 4: (a) The phase diagram obtained from the self- 
consistent numerical iterations of the gap equations for \m = 
0.52, As = 0.30. The tetracritical point corresponds to the 
point C in the Fig. [T] The green shade in fig(a) represents 
the coexisting region of IC-SDW and SC. The dashed line 
branching off from the T^j represents the boundary between 
C and IC-SDW. In figures (b) and (c), the order parameters 
along the dotted fines of figure (a) are plotted. The red fine 
is A and the blue line is m as a function of temperature at 
5 = 0.18 and as a function of doping T = 0.019, respectively. 
All transitions below T^ are continuous. Note in figure (b) 
that m decreases as T is lowered below the coexisting region 
because SC begins to gap out the Fermi surface partially. 



ary. This makes the phase boundaries more complicated 
than previously reported as presented in Fig. 1. 

Finally, we remark that the shapes of the electron 
Fermi surfaces take quite different forms for different 
class of the pnictides. The effect of the FS shapes 
on the phase coexistence was studied in RefjiS,. Also, 
the effect of pairing symmetry on the phase separa- 
tion/coexistence will be interesting. We are currently ap- 
plying the present method to understand how the phase 
coexistence/separation phenomenon depends on the pair- 
ing symmetry. 
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incommensurate SDW cases. The energy of Eq. p5l 
above is generahzed for the non-zero incommensurabil 

ity (q ^ 0) to 



Appendix A: Incommensurate SDW 

In the appendix, we cohect the technical details of the 
Ginzburg-Landau expansion of the free energy for the 
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yrn^ + {£, — q cos 6) ± {6 + q cos ( 




and the SDW order parameter m and superconducting order parameter A to 
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The 4th order coefficients, /3a, Pm, and /?', of the Ginzburg-Landau expansion of the free energy around the 
multicritical point of Eq. p3)) may be obtained by the derivative of the free energy of Eq. (fT6|) with respect to the 
order parameters. They are given by 
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(A6) 



r 



The 4th order coefficients are then used to compute the and (|35p . The /3 determines the phase coexistence or 
anisotropy coefficient /3 using the equations ((^ , ([50)1 , separation by the condition of Eq. ((551) . 
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